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O ' ABSTRACT 

| We use our most recent training set for the Rico code to estimate the impact of recom- 

bination uncertainties on the posterior probability distributions which will be obtained 

■ from future CMB experiments, and in particular the Planck satellite. Using a Monte 

Carlo Markov Chain analysis to sample the posterior distribution of the cosmological 

Q • parameters, we find that Planck will have biases of —0.7, —0.3 and —0.4 sigmas for 

pj ' us, flbh 2 and log(10 10 As), respectively, in the minimal six parameter ACDM model, if 

c/2 , the description of the recombination history given by RiCO is not used. The remaining 

parameters (e.g. r or Qa m h 2 ) are not significantly affected. We also show, that the 
cosmology dependence of the corrections to the recombination history modeled with 
, RiCO has a negligible impact on the posterior distributions obtained for the case of the 

J> ■ Planck satellite. In practice, this implies that the inclusion of additional corrections 

£f~) [ to existing recombination codes can be achieved using simple cosmology-independent 

OO . 'fudge functions'. 

CO ■ Finally, we also investigated the impact of some recent improvements in the treat- 

ment of hydrogen recombination which are still not included in the current version of 
. our training set for RiCO, by assuming that the cosmology dependence of those correc- 

t-H ' tions can be neglected. In summary, with our current understanding of the complete 

OS , recombination process, the expected biases in the cosmological parameters inferred 

■ from Planck might be as large as —2.3, —1.7 and —1 sigmas for ng, Vt^h 2 and 
| log(10 10 As) respectively, if all those corrections are not taken into account. We note 

• i-H . that although the list of physical processes that could be of importance for Planck 

1 seems to be nearly complete, still some effort has to be put in the validation of the 

J-j \ results obtained by the different groups. 

. . . i The new RiCO training set as well as the fudge functions used for this paper are 

publicly availabe in the RiCO-webpage. 
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1 INTRODUCTION 

The cosmic microwave background (CMB) is nowadays an 
essential tool of theoretical and observational cosmology. Re- 
cent advances in the observations of the CMB angular fluctu- 
ations in temperature and polarization (e.g. iHinshaw et al.l 
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f E-mail: jchlubaOcita.utoronto.ca 



l2009l ; lKomatsu et alj|2009h provide a detailed description of 
the global properties of the Universe, and the cosmological 
parameters are currently known with accuracies of the or- 
der of few percent in many cases. The experimental prospec t 
for the Planck satellite (|The Planck Collaboration! 120061 ) , 
which was launched on May 14th 2009, is to achieve the 
most detailed picture of the CMB anisotropies down to an- 
gular scales of I ~ 2500 in temperature and i ~ 1500 in 
polarization. This data will achieve sub-percent precision in 
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many cosmological parameters. However, those high accura- 
cies will rely on a highly precise description of the theoret- 
ical predictions for the different cosmological models. Cur- 
rently, it is widely recognised that the major limiting factor 
in the accuracy of angular power spectrum calculations is 
the uncertainty in the ionization h istory of the Universe (see 
iHu et alJ[l99a ISeliak et alj|2003h . 

This has motivated several groups to re-examine the 
prob l em of cosmol ogical recombination l)Zeldovich et al.l 
1 19681 : |Peebledfl96sh . taking into account detailed correc- 
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corrections individually leads to changes in the ionization 
history at the level of > 0.1%, in such a way that the cor- 
responding overall uncertainty in the CMB angular power 
spectra exce eds the benchmark of ±3/£ at large i (for more 
details, see iFendt et all 120091 . hereafter FCRW09), thus 
biasing any parameter constraints inferred by experiments 
like Planck, which will be cosmic variance limited up to 
very high multipoles. 

The standard description of the recombin ation process 
is pro vided by the widely used Recfast code ijSeaeer et al.l 
119991 ). which uses effective three-level atoms, both for hydro- 
gen and helium, with the inclusion of a conveniently chosen 
fudge factor which artificially modifies the dynamics of the 
process to reproduce the resu lts of a multilevel recombina- 
tion code (jSeaeer et al.ll2000r i. 

The simultaneous evaluation of all the new effects dis- 
cussed above make the numerical computations very time- 
consuming, as they currently require the solution of the full 
multilevel recombination code. Moreover, some of the key 
ingredients in the accurate evaluation of the recombination 
history (e.g. the problem of radiative transfer in hydrogen 
and the proper inclusion of two-photon processes) are solved 
using computationally demanding approaches, alth ough in 
some cases semi-analytical approximations (see e.g. iHiratal 
120081 ') might open the possibility of a more efficient evalua- 
tion in the future. 

In order to have an accurate and fast representation of 
the cosmological recombination history as a function of the 
cosmological parameters, two possible approaches have been 
considered in the literature. The first one consists of the in- 
clusion of addition al fudge factors to mimic the new physics, 
as recently done in IWong et al] (|2008h (see Recfast vl.4.2), 
where they include an additional fudge factor to modify the 
dynamics of helium recombination. The second approach is 
the so-called RlCO code (FCRW09), which provides an ac- 
curate representation of the recombinati on history by us- 
ing a regression scheme based on Pico ()Fendt fc WandeTtl 
l2007al lbh. The RlCO code smoothly interpolates the X e (z;p) 
function on a set of pre-computed recombination histories 
for different cosmologies, where 2 is the redshift and p rep- 
resents the set of cosmological parameters. 



In this paper, we present the results for parameter es- 
timations using RlCO with the most recent training set pre- 
sented in FCRW09. This permits us to accurately account 
for the full cosmological dependence of the corrections to the 
recombination history that were included in the multi-level 
recombination code which was used for the training of RlCO 
(see Sect. l2.l1 for more details). With this tool, we have eval- 
uated the impact of the corrections on the posterior proba- 
bility distributions that are expected to be obtained for the 
Planck satellite, by performing a complete Monte Carlo 
Markov Chain analysis. The study of these posteriors have 
shown that the impact of the cosmology dependence is not 
very relevant for those processes included into the current 
RlCO training set. Therefore, by assuming that the cosmol- 
ogy dependence of the correction in general can be neglected, 
we have also investigated the impact of recent improvements 
in the treatment of hydrogen recombination (see Sect. 12.2")) . 
The basic conclusion is that, with our current understand- 
ing of the recombination process, the expected biases in the 
cosmological parameters inferred from Planck might be as 
large as 1.5-2.5 sigmas for some parameters as the baryon 
density or the primordial spectral index of scalar fluctua- 
tions, if all these corrections to the recombination history 
are neglected. 

The paper is organized as follows. Sect. [2] describes the 
current training set for RlCO, and provides an updated list 
of physical processes during recombination which were not 
included in FCRW09. Sect. [3] presents the impact of the 
recombination uncertainties on cosmological parameter esti- 
mation, focusing on the case of Planck satellite. Sect.[4]fur- 
ther extends this study to account for the remaining recom- 
bination uncertainties described in Sect. [5] Sect. O presents 
the analysis of present-day CMB experiments, for which the 
effect is shown to be negligible. Finally, the discussion and 
conclusions are presented in sections [S] and [7] respectively. 



2 UPDATED LIST OF PHYSICAL PROCESSES 
DURING RECOMBINATION 

In this Section we provide an updated overview on the im- 
portant physical processes during cosmological recombina- 
tion which have been discussed in the literature so far. We 
start with a short summary of those processes which are 
already included into the current training set (FCRW09) of 
RlCO fSect. l2~Tj) . The corresponding correction to the ioniza- 
tion history close to the maximum of the Thomson visibility 
function is shown in Fig. [1] 

We then explain the main recent advances in connection 
with the radiative transfer calculations during hydrogen re- 
combination (Sect. I2.2|l . which lead to another important 
correction to the cosmological ionization history (see Fig. [TJ 
that is not yet included into the current training set of RlCO. 
However, as we explain below (Sect. 3.4) it is possible to take 
these corrections into account (Sect. 4), provided that their 
cosmology dependence is negligible. Our computations show 
that this may be a valid approximation (Sect. 3.4). 

We end this section mentioning a few processes that 
have been recently addressed but seem to be of minor impor- 
tance in connection with parameter estimations for Planck. 
Overall it seems that the list of processes that could be of 
importance in connection with Planck is nearly completed. 
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Figure 1. Recombination uncertainties with respect to the stan- 
dard Recfast vl.4.2 code. For the fiducial cosmological model, 
it is shown the correction to the recombination history which is 
incorporated in the new RlCO training set (d otted line); the cor- 
rection due to Lya radiative transfer effects l lChluba fc bunyacv 
2009ajQI3) (dot-dashed); the effect of Raman scattering ijHirata 
20081) (dashed); and the combination of all previous effects (solid 



Figure 2. Same as figure [T] but for the angular power spectrum 
(APS). For the fiducial cosmological model, it is shown the cor- 
rection to the APS due to the recombination history which is 
incorporated in the new RlCO training set (gray lines), and the 
correction due to the additional inclu sion of Lya radiative trans- 
fer effects JChluba fc bunyac vl l2009allclle] ) and Raman scattering 
i]Hiratall2008l) (dark lines). 



However, still some effort has to go in cross-validation of the 
results obtained by different groups. 



2.1 The current training set for RlCO 

As demonstrated in FCRW09, RlCO can be used to repre- 
sent the recombination history of the Universe, accurately 
capturing the full cosmology dependence and physical model 
of the multilevel recombination code that was used in the 
computations of the RlCO training set. 

For the current RlCO training set we ran our full re- 
combination code using a 75-shell model for the hydrogen 
atom. The physical processes which are included during hy- 
drogen recombination are described in detail in FCRW09: 
the in duced 2s- Is two-photon decay ()Chluba fe Sunvaevl 
l2006bl ); the feed back of Lyman a photons on the ls-2s 
absorption rate ()Kholupenko fc Ivanchikl l2006t ) ; the non- 
equilib rium populations in the angular momentum sub- 
states (|Rubino-Martm et al.ll2006l ; IChluba et alj|2007l); and 
the effect of Lyman series feedback ( Chluba fc Sunvaevl 
120071 ). For helium recombination we took into ac- 



count: the spin- forbidden H e I 2 3 Pi — l 1 So transition 
(jDubrovich fc Grachevl 120051 ): and the acceleration of he- 
lium r ecombination by neutral hydrogen |Switzer fc Hiratal 
l2008al ). Furthermore, we also updated our physical constants 
according to the NlST database!, including the new value of 
the gr avitational constant and the helium to hydrogen mass 
ratio ijWong fc Scottll2007l ). 

A more detailed description of the physical processes 
that were taken into account in the current RlCO training 
set can be found in FCRW09. This RlCO training set is now 
publicly available at http : //cosmos . astro .uiuc . edu/rico. 



2.2 Updated radiative transfer calculations during 
hydrogen recombination 

As al ready pointed out earlier (e.g. IChluba fc Sunvaevl 
l2008bl ) in particular a detailed treatment of the hydrogen 
Lyman a radiative transfer problem including two-photon 
corrections is expected to lead to an important additional 
modification during hydrogen recombination. An overview 
of the relevant physical aspects in connection with this prob- 
lem was already given in Sect. 2.3.2 and 2.3.3 of FCRW09. 
However, at that time the problem was still not solved at 
full depth, but recently several important steps were taken, 
which here we n ow want to discuss briefl y (for additional 
overview see also ISunvaev fc Chlubal 120091 ). 



2.2.1 Partial frequency redistribution due to line 
scattering 

Recently, the effect of partial frequency redistribution on 
the Lyman a escape rate and the ionization history during 
hy drogen recombination was i ndep endently studied in deta il 
bv lChluba fc Sunvaevl (l2009al) andlfflrata fc Forbes! < |200Sh . 
As shown by Chluba fc Sunvaevl (!2009aT l. the atomic recoil 
effect leads to the dominant contribution to the associated 
correction in the ionization history, and the result for this 
process alone see ms to be in good agreement w ith the one 
obtained earlier bv lGrachev fc Dubrovichl l|2008h . leading to 
AN C /AN C ~ -1 - 2% at z ~ 900. Also the computations by 
iHirata fc Forbes! (120091') seem to support this conclusion. 

However. IChluba fc Sunvaevl (|2009al ) also included the 
effect of Doppler broadening and D oppler boosting which 
was n eglected in the analysis of iGrachev fc Dubrovichl 
l|2008l ). Doppler boosting acts in the opposite direction as 
atomic recoil and therefore decelerates recombination, while 
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2 They used a Fokker-Planck approximation (e.g. lRvbickill2006t ) 
for the frequency redistribution function. 
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the effect of Doppler broadening can lead to both an in- 
crease in the photons escape probability or a decrease, 
depending on the initial frequency of the photons (see 
IChluba fc Sunvaevl [2009al . for more detailed explanation). 
The overall correction to the recombination history due to 
line scatterin g amounts to AN B / AN C ~ — 0.6% at z ~ 900. 
The results of IChluba fc Sunvaevl (1200981) seem to be rather 
similar to those of Hirata fc Forbe si (|2009l ). however a fi- 
nal comparison will become necessary in order to reach full 
agreement on the final correction. In th e computations pre- 
sented below we will use the results of IChluba fc Sunvaevl 
l|2009ah . 



2.2.2 Two-photon transitions from higher levels 

Initially, the problem of two-photon transitions from highly 
excited levels in hydrogen and helium was proposed by 
iDubrovich fc Grachev! (|2005h . However, for hydrogen recom- 
bination onl y very recently this problem has been solved con - 
vincingly bv lHiratal (|2008l ) and lChluba fc Sunvaevl (|2009ch . 
using two inde pendent, conceptually differ ent approaches. 
Also until now, IChluba fc Sunvaevl (|2009d ) took the main 
contribution coming from the 3d- Is and 3s- Is two -photon 
profile corrections into account, while IHiratal (|2008h also in- 
cluded the ns-ls and nd-ls two-photon profile corrections 

for larger n. 

In the analysis of lChluba fc Sunvaevl (|2009(J ). three in- 
dependent sources for the corrections in connection with the 
two-photon picture were identified (we will discuss the other 
two processes in Sect. I2.2~4l and l2.2,3p . As they explain, the 
total modification coming from purely quantum mechanical 
aspects of the problem (i.e. corrections due to deviations of 
the lin e profiles from the normal L orenzian shape, as pointed 
out bv lChluba fc Sunvaevl (|2008bl l) leads to a change in the 
free electron number of AN e /AN c ~ -0.4% at 2 ~ 1100. 
It also seems clear that the remaining small difference (at 
th e level of AN e /N e ~ 0.1% - 0.2%) between the results 
of lHiratal l|200g| ) and lChluba fc Sunvaevl (|2009cl ) for the to- 
tal correction rel ated to the two-phot o n deca ys from excited 
states is because IChluba fc Sunvaevl (|2009cl ) only included 
the full two-photon profiles of t he 3s-ls and 3d-ls channels . 
Below we will use the results of lChluba fc Sunvaevl (|2009cT l 
for this process. 



2.2.3 Time- dependent aspects in the emission and 
absorption of Lyman a photon 

One of the key ingredients for the derivation of the escape 
probability in t he Lyman a resonance using the Sobolev 
approximation (|Sobolevl li"96bT l is the quasi- stationarity of 
the line transfer problem . However, as shown recently 
(|Chluba fc Sunvaevl l2009el lch at the percent-level this as- 
sumption is not justified during the recombination of hy- 
drogen, since (i) the ionization degree, expansion rate of 
the Universe and Lyman a death probability change over a 
characteristic time Az/z ~ 10%, and (ii) because a signifi- 
cant contribution to the total escape probability is coming 
from photons emitted in the distant wings (comparable to 
10 2 - 10 3 Doppler width) of the Lyman a resonance. There- 
fore, one has to include time-dependent aspects in the emis- 
sion and absorption process into the line transfer problem, 



leading to a delay of recombination by ANe/Ne ~ +1.2% at 
z ~ 100 0. Below we will use the results of lChluba fc Sunvaevl 
l|2009ch . 



2.2-4 Thermodynamic asymmetry in the Lyman a 
emission and absorption profile 

As explained bv lChluba fc Sunvaevl (|2009d ). the largest cor- 
rection related to the two-photon formulation of the Lyman 
a transfer problem is due to the frequency-dependent asym- 
metry between the emission and absorption profile around 
the Lyman a resonance. This asymmetry is given by a 
thermodynamic correction factor, which has an exponen- 
tial dependence on the detuning from the line center, i.e. 
/„ oc exp[h(u — u a )/kT-y\, where v a is the transition fre- 
quency for the Lyman a resonance. Usually this factor can 
be neglected, since for most astrophysical problems the main 
contribution to the number of photons is coming from within 
a few Doppler width of the line center, where the thermo- 
dynamic factor indeed is very close to unity. However, in 
the Lyman a escape problem during hydrogen recombina- 
tion contributio ns from the very distant da mping wings are 
also important (|Chluba fc Sunvaev 2009c. c), so that there 

f v 7^ 1 has to be incl uded. 

As explained in (|Chluba fc Sunvaevl l2008bh . the ther- 
modynamic factor also can be obtained in the classical 
picture, using the detailed balance principle. However, in 
the two-photon picture this factor has a natural explana- 
tion in connection with the absorption of photons from the 
CMB blackbody ambient radiatio n field (|Chluba fc Sunvaevl 
l2009d : ISunvaev fc Chlubal l2009h . This process leads to a 
~ 10% increase in the Lyman a escape probability, and 
hence accelerates hydrogen rec ombination. For the correc- 
tion to the ionization history, IChluba fc Sunvaevl (|2009ch 
obtained AN c /A N e ~ -1.9% at z ~ 1100. Note also that 
in the analysis of lChluba fc Sunvae v ( 2009c) the thermody- 
namic correction factor was included for all ns-ls and nd-ls 
channels with 3 ^ n ^ 10 . Below we will use the results of 
IChluba fc Sunvaevl (|2009d ) for this process. 



2.2.5 Raman scatterings 

iHiratal l|2008l ) also studied the effect of Raman scatterings 
on the recombination dynamics, leading to an additional 
delay of hydrogen recombination by AN e /N c ~ 0.9% at z ~ 
900. Here in particular the correction due to 2s-ls Raman 
scatterings is important. Again it is expected that a large 
part of this correction can be attributed to time-dependent 
aspects and the correct formulation using detailed balance, 
and we are currently investigating this process in detail. In 
the computat ions presented below we will use the results of 
IHiratal (|2008l ) for the effect of Raman scatterings. 



2.3 Additional processes 

There are a few more processes that here we only want men- 
tion very briefly (although with this the list is not meant to 
be absolutely final or complete), and which we did not ac- 
count for in the computations presented here. However, it is 
expected that their contribution will not be very important. 
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2. 3. 1 Effect of electron scattering 

The effect of electron scattering during hydrogen recombi- 
nation was also recently investigated by IChluba fc Sunvaevl 
(|2009ch using a Fokker-Planck approach. This approxima- 
tion for the frequency redistribution function may not be 
sufficient towards the end of hydrogen recombination, but 
in the overall correction to the ionization history was very 
small close the maximum of the Thomson visibility func- 
tion, so that no big difference are expected when more ac- 
curately using a scatt ering Kernel-approach. Very recently 
lHaimoud et al.l (|2009T ) showed that this statement indeed 
seems to be correct. 



2.3.2 Feedback of helium photons 

Very recently IChluba fc Sunvaevl (|2009bf ) investigated the 
feedback problem of helium photons including the processes 
of 7(Hei) — > He i, 7(Hei) — > Hi, 7(Heii) — » Hei and 
7(Heii) — > Hi feedback. They found that only 7(Hei) — > 
Hei feedback leads to some small correction (AN C /N C ~ 
+0.17% at z ~ 2300) in the ionization history, while all 
the other helium feedback induced corrections are negligi- 
ble. This is because the 7(Hei) — > Hi, 7(Heii) — > Hei 
and 7(Heii) — > Hi feedback processes all occur in the pre- 
recombinational epochs of the considered species, where the 
populations of the levels are practically in full equilibrium 
with the free electrons and ions. 

T he 7(Hei) — > Hei feedbac k process was already stud- 
ied by Switzcr & Hirata (2008a|), but the result obtained by 
IChluba fc Sunvaevl (|2009bf l seems to be smaller. However, 
it is clear that any discrepancy in the helium recombination 
history at the 0.1% - 0.2% level will not be very important 
for the analysis of future CMB data. 

We would also like to mention, that although the 
7(Hei) — > Hi, 7(Heii) — » Hei and 7(Heii) — > Hi feedback 
processes do not affect the ionization history, they do intro- 
duce interesting changes in the recombinational radiation, 
increasing t he total contribution of ph otons from helium by 
40% - 70% (jChluba fc Sunvaevll2009bh . 



2.3.3 Other small correction 

Recently, several additional processes during hydrogen re- 
combination were discussed. These include: the overlap of 
Lyman ser ies resonances caused by the thermal motion of 
l2009h ; 



the atoms 



Hai moud et al 
— T — 



hydroge n (IGrin et al.l 1200911: hydrogen deuterium recombi- 
nation (|Fung fc Chlubal 120091 : iKholupenkol 120091); and 3s- 
2s and 3d-2s two-photon transitions ( Chluba fc Sunvaevl 



2008b). All these processes seem to affect the ionization his- 
tory of the Universe at a level (well) below 0.1%. 

Furthermore, one should include the small re-absorption 
of photons from the 2s- Is two- photon continuum close to 
the Lyman a resonance, where our estimates show that this 
leads to another AN e /N e ~ 0.1% - 0.2% correction. 

2.4 Overall correction 

Figure [l] shows the current best-estimate of the remain- 
ing overall correction to the recombination history, while 



Fig- [U translates these corrections into changes of the angu- 
lar power spectrum. To descri be the Ly-q transfer effect s, 
we adopt the curve presented in IChluba fc Sunvaevl (|2009ah . 
which includes the results of t he processes investigated b y 
I Chluba fc Sunvaevl (|2009<J ) and lChluba fc Sunvaevl (|2009el ). 
We also include the effect of Raman scatter ing on hydrogen 
recombination as described in iHiratal (|2008l ). 

It seems that the remaining uncertainty due to pro- 
cesses that were not taken into account here can still exceed 
the 0.1% level, but likely will not lead to any significant 
addition anymore. However, it is clear that a final rigorous 
comparison of the total result from different independent 
groups will become necessary to assure that the accuracy 
required for the analysis of Planck data will be reached. 
Such detailed code comparison is currently under discussion. 



3 IMPACT OF RECOMBINATION 

UNCERTAINTIES ON COSMOLOGICAL 
PARAMETER ESTIMATION 

Here, the RlCO code together with the latest training set de- 
scribed above is used to evaluate the impact of the recom- 
bination uncertainties in the cosmological parameter con- 
straints inferred from CMB experiments. 

All our anal yses use the soft ware packages CAM^ fl 
(|Lewis et alJ^OOfj ) and CosmoMcQ (jLewis fc Bridlell2002r i. 
Camb is used to calculate the linear-theory CMB angu- 
lar power spectrum, and has been modified here to include 
the recombination history as described by RlCO. CosmoMC 
uses a Monte Carlo Markov Chain (MCMC) method to sam- 
ple the posterior distribution of the cosmological parameters 
from a given likelihood function which describes the exper- 
imental constraints on the CMB angular power spectrum. 

The default parametrization which is included inside 
CosmoMC exploits some of the intrinsic degeneracies in 
the C MB angular power spectrum (see e.g. iKosowskv et al.l 
120021 ) . and uses as basic parameters the following subset: 



Patd = {Qbh , n dm /i ,0,t, n s ,log(10 A s )}. 



(1) 



Here, fit,/i 2 and f2d m /i 2 are the (physical) baryon and dark 
matter densities respectively, where h stands for the Hub- 
ble constant in units of 100 km s _1 Mpc -1 ; 8 is the acous- 
tic horizon angular scale; r is the Thomson optical depth 
to reionization; and ns and As are the spectral index and 
the amplitude of the primordial (adiabatic) scalar curvature 
perturbation power spectrum at a certain scale fco. For the 
computations in this paper, we will use fco = 0.05 Mpc -1 . 

It is important to note that the above parametrization 
p s td ma kes use of an approxi mate formula for the sound 
horizon (|Hu fc Sugiv ama 199^), which used in its derivation 
some knowledge on the recombination history provided by 
Recfast. Therefore, the use of this parameter is not appro- 
priate if we are changing the recombination history, as this 
might introduce artificial biases. For this reason, in this pa- 
per we have modified CosmoMC to use a new parametriza- 
tion pkew, which is defined as 

Pnew = {O b /l 2 , fldmh 2 , H , T, US , log( 1Q 10 A S ) } . (2) 



3 http://camb.info/ 

4 http:/ /cosmologist.info/cosmomc/ 
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where we have replaced 9 by Ho as a basic parameter. In 
this way, we still exploit some of the well-known degenera- 
cies (e.g. that between r and As), but we eliminate the 
possible uncertainties which may arise from the use of 6, 
at the expense of decreasing slightly the speed at which 
the chain converges. As we show below, both the shape of 
the posteriors and the confidence levels for the rest of the 
parameters is practically unaffected by this modification in 
the parametrization, despite of the fact that the CosmoMC 
code now assumes a flat prior on Ho instead of a flat prior 
on 9. A modified version of the params.f90 subroutine in- 
side CosmoMC which uses this new parametrization is also 
available in the RlCO webpage. 

Finally, concerning the convergence of the chains, all 
computations throughout this paper were obtained using 
at least five independent chains; th ose chains have been 
run until the lGelman fc Rubinl l| 19921 ) convergence criterion 
R — 1 yields a value smaller than 0.005 for the minimal (6- 
parameter) case, and values smaller than 0.02 — 0.2 for those 
cases with a larger number of parameters. 



3.1 Impact on parameter estimates for Planck 
alone 

For the case of Planck satellite, the mock data is prepared 
as follows. We assume a Gaussian symmetric beam, and the 
noise is taken to be uniform across the sky. For definiteness, 
we have adopted the nominal values of the beam and pixel 
noise whic h correspond to the 143 GHz Pla nck band, as de- 
scribed in lThe Planck Collaboration! (|2006h . Thus, we have 

#beam = 7.1', W^ 1 = <7 2 oisc r2 bcam = 1-53 X 10~ 4 /iK 2 and 

Wp 1 = 5.59 x 10~ 4 /iK 2 , where wt and Wp indicate the in- 
tensity and polarization sensitivities, respectively. Note that 
using a single frequency channel is implicitly assuming that 
the remaining Planck frequencies have been used to fully 
remove the foreground contamination from this reference 
143 GHz channel. For a complete discussion on which com- 
bination of Planck channels is more appropriate in t erms 
of the parameter constraints, see lColombo et all (|2009l ). 

The fiducial cosmological model used for these 
computations corres ponds to the WMAP5 cosmology 
(jHinshaw et al.ll2009f ). and has parameters Q,bh 2 = 0.02273, 
ftdm/i 2 = 0.1099, h = 71.9, t = 0.087, n s = 0.963 and an 
amplitude 2.41 x 10~ 9 at k = 0.002 Mpc -1 (or equivalently, 
A s = 2.14 x 10~ 9 at fc = 0.05 Mpc -1 , or log(10 10 A s ) = 
3.063). The mock data is then produced using the cosmolog- 
ical recombination history as computed with our complete 
multi-level recombination code for this particular cosmol- 
ogy (using n = 75 shells to model the hydrogen atom). The 
shape of the likelihood function adopted here corresponds 
to the exact full-sky Gaussian likelihood function given in 
iLewis! (|2005h . which is implemented in CosmoMC using the 
all_l_exact format for the dat£0. We note that for the 
mock Planck data we use not a simulation but the actual 
(fiducial) angular power spectrum. 

We run CosmoMC for the case of Planck mock data 
alone and the minimal model with six free parameters. For 



5 See also http://cosmocofFcc.info/vicwtopic. php?t=231 for more 
details. 



each set of runs, we consider two cases for the recombina- 
tion history, one is the current version of Recfast (vl.4.2), 
and the second one is the RlCO code with our latest training 
set. The main result are summarized in table [T] and the cor- 
responding posterior distributions are shown in Fig. [3] We 
note that the sizes of the error bars on these parameters are 
similar to those obtained for Planck by other authors (e.g 
iBond et al1l2004UColombo et al.ll2009r i. For this six param- 
eter case, the largest biases do appear in ns (—0.7 sigmas), 
n b h 2 (-0.3 sigmas) and log(10 10 A s ) (-0.4 sigmas). 

The sign of the correction in ns, which is the param- 
eter having the largest bias, can be understood as follows. 
The physical effects to the recombination history included 
in RlCO produce a slight delay of the recombination around 
the peak of the visibility function, i.e. an excess of electrons 
with respect to the standard computation (see Fig. [TJ. This 
in turn produces a slightly larger Thomson optical depth, 
which increases the damping of the anisotropies at high mul- 
tipoles. In order to compensate this excess of damping, the 
analysis which uses the standard Recfast code gives a lower 
value of ns. 

The biases on the other three parameters, i.e. r (which 
is well-constrained by large scale polarization measure- 
ments), Qdmh 2 and Ho, are negligible (i.e. less than 0.1 sig- 
mas). However, if the analysis is repeated with the standard 
parametrization using 9 as basic parameter, then we would 
find a +1.8 sigma bias in 9, while for the rest of the parame- 
ters the posteriors remain the same as before. For illustration 
of these facts, Fig. fl] shows a comparison of the posteriors 
for ns obtained with the two different parametrizations, as 
well as the posteriors obtained for 9 if using the standard 
parametrization. 

Finally, we have checked that the inclusion of lensed 
angular power spectra on the complete procedure modifies 
neither the shape of the posteriors nor the biases. Therefore, 
for the rest of the paper we perform the computations in the 
case without lensing. This decreases the computational time 
by a significant factor. 

3.2 Importance of the cosmology dependence of 
the correction 

One of the questions that can be explored with RlCO and 
the new training-set is the importance of the cosmology de- 
pendence of the corrections to the recombination history. 
In order to obtain a simplified description of the recombi- 
nation history, it is important to evaluate if the cosmology 
dependence of the corrections plays a significant role in de- 
termining the final shape of the posteriors. To explore this 
issue, we have modified the standard Recfast code by in- 
troducing the following (cosmology independent) correction: 

x" cw (z; {cosmology}) = xf ECFAST (z; {cosmology})/^) 

(3) 

where this function f(z) is computed as 

/(*) = ! +^2 (4) 

for a certain fiducial cosmological model. By introducing 
this modification inside CosmoMC, we have compared the 
posterior distributions obtained in this case with those ob- 
tained using the full RlCO code with the new training set. 
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Figure 3. Impact on parameter estimates for Planck satellite, in the six-parameter case. Dark solid line shows the 1-D and 2-D posterior 
distributions which are obtained when the RlCO code is used to describe the recombination history. Light dashed solid line shows the 
same posteriors obtained with the Recfast code. The biases are evident in rig, Q h h 2 and log(10 10 A s ). 



Table 1. Impact of recombination uncertainties on the confidence limits of the cosmological parameters for the Planck satellite. Confidence 
intervals are derived as the 0.16, 0.5 and 0.84 points of the cumulative probability distribution function, in such a way that our parameter 
estimate is the median of the marginalised posterior probability distribution function, and the confidence interval encompasses 68 per 
cent of the probability. The biases correspond to the case of using the Rico training set. 





With Rico 


With Recfast 
(vl.4.2) 


Bias 
(in sigmas) 


Fiducial 
model 


n h h 2 (xio 2 ) 
n dm h 2 

Ho 

T 

n s 

log(10 10 A s ) 


2.273 ±0.014 
n inq8+ - 0013 

o.ioys„ 0012 

71 9+ ' 6 
0.087 ± 0.006 

a 6 q 2 + 0038 
u - 3uo -0.0034 
, 06 o+0.009 
°- UUJ -0.008 


2.269 ±0.014 

1099+ ' 0012 
u.j.u»»_ 0013 

72 0+ - 7 

0.087 ±0.006 

9606+ 0035 

U.9DUD_ 0037 

3.059 ± 0.009 


-0.31 
0.06 
0.11 
-0.04 
-0.74 
-0.42 


2.273 
0.1099 
71.9 
0.087 
0.963 
3.063 



The results of this comparison are shown in Figure [5] The 
fact that the differences in the posterior are negligible indi- 
cates that there is no significant cosmological dependence in 
the corrections to the ionization history included in the RlCO 
training set. Therefore, and for the case of the Planck satel- 
lite, one can in principle use cosmology independent 'fudge 
functions', as the one presented in Eq. to accommodate 
additional corrections to the recombination history. 

As a final check in this section, we have explored the 
sensitivity of the fudge function f(z) to the cosmological 
model which is used as fiducial model for the full compu- 
tation of the recombination history. Taking as a reasonable 
range of variation the two-sigma confidence interval which 



is ob tained of the analysis of WMAP5 data (|Dunklev et al.l 
l2009t ). we have compared the f{z) functions obtained from 
cosmological models which differ two-sigmas with respect to 
the actual fiducial model which is used in this paper. The 
result is that the changes in f(z) with respect to the solid 
curve presented in Fig. [T]are below one percent, thus giving 
a negligible correction to the main effect included in f(z). 
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Figure 4. Top: Posterior distributions for ng obtained with 
two different parametrizations. The solid line corresponds to the 
marginalised 1-d posterior distribution obtained using RlCO while 
the dashed line corresponds to the case of using Recfast, both 
with the new parametrization. The dotted lines show in each 
case the posteriors obtained when the standard parametrization 
is used (i.e. 8 instead of Hq as fundamental parameter). Bottom: 
Bias recovered on 6 parameter if the standard parametrization is 
used. This case corresponds to the same case as Fig. [3] i.e. six 
parameter case. 



4 ESTIMATING THE CORRECTIONS FROM 
THE REMAINING RECOMBINATION 
UNCERTAINTIES 

Although a full recombination code which includes all the 
physical effects discussed in Sect.[2]is still not available, there 
is a good agreement in the community about the list of rele- 
vant physical processes that have to be included. Moreover, 
all those effects have been already discussed in the literature 
by at least one group, so we have estimates of the final im- 
pact of these corrections on the recombination history (see 

Fig. ID. 

Based on these estimates, we have quantified the im- 
pact of future corrections to the recombination history using 
the approximation described in Eq. ((3]), where f(z) is taken 
from Fig. [T] as described in Eq. Q. Using this function, 
we have obtained the posterior distributions for the same 
case discussed above (nominal Planck satellite sensitivities 
and a six parameter analysis). The basic results are shown 
in Table [5] and Fig. [H] The biases in the different param- 
eters increase very significantly, as one would expect from 
the inspection of Fig. [T] specially for ns (—2.3 sigmas), Qhh 2 
( — 1.65 sigmas) and log(10 10 As) (—1 sigmas). Therefore, as 




0.955 0.960 



0.965 0.970 0.975 




0.0224 0.0226 0.0228 0.0230 




Figure 5. Comparison of the 1-D posterior distributions obtained 
with the full RlCO code and the new training set (solid lines), with 
those obtained with the simplified description given in Eq. (|3jl 
(dashed lines). For this case, the fiducial model which has been 
used to compute f(z) is taken to be the same as the fiducial 
cosmology. In all four panels, the vertical dotted line shows the 
value of the fiducial model. 



pointed out by seve ral authors (see e.g. IChluba fc Sunvaevl 
l2008bl ; lHiratall2008l ) . the detailed treatment of the hydrogen 
Lyman a radiative transfer problem constitutes the most 
significant correction to our present understanding of the 
recombinational problem. If this effect is not taken into ac- 
count when analysing Planck data, the final constraints 
could be significantly biased. 

4.1 Extended cosmological models 

In this subsection we describe to what extend the full set of 
corrections to the recombination history may affect the cos- 
mological constraints on some extensions to the (minimal) 
six-parameter model which was used in this work. Through- 
out this subsection, we compare the complete recombination 
history (which includes the additional corrections shown in 
Fig. [TJ) with the constraints that would be inferred using 
Recfast vl.4.2. 

4-1.1 Tensor perturbations 

We first consider the case of including tensor perturbations 
in addition to the previous model. For these computations, 
we consider an 8-parameter model, by including the spectral 
index of the primordial tensor perturbation (tit) and the 
tensor-to-scalar ratio r, in addition to the parameters in 
equation As pivot scale, we are using fco = 0.05 Mpc -1 . 

Fig. [7] shows the impact of recombination uncertainties 
on the r-ns plane. The constraints on r are determined by 
large-scale information, and therefore the modifications of 
the recombination history do not bias this parameter. How- 
ever, due to the important bias on ns, the 2-D contours on 
this plane are significantly shifted. 

4-1.2 Scale dependence of spectral index 

The running of the spectral index is a possible extension to 
the simple ACDM model which is under debate in the liter- 
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Table 2. Estimate of the impact of remaining recombination uncertainties on the confidence limits of the cosmological parameters for 
the Planck satellite. Confidence intervals are derived as in tabled so the confidence interval encompasses 68 per cent of the probability. 
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Figure 6. Posterior distributions obtained with the remaining recombination uncertainties (see text for details). Solid line corresponds 
to the posteriors using the correct description of the recombination history, while the dashed lines represent the case in which the current 
description (Recfast vl.4.2) is used. 



ature in light of WMAP observations (see e.g. ISpereel et al.l 
,20071 . which gives n run = dn s /d\nk = -0.055±0.030). Fig.[8] 
illustrate the impact of the recombination uncertainties on 
the determination of the running of the spectral index. For 
this computation, we considered a 7-parameter model by 
adding n run = dns/dlnk, which again is computed at the 
same pivot scale of ko — 0.05 Mpc -1 . When studying the 
one-dimensional posterior distributions for all parameters, 
we find that the inclusion of n run does not affect the shape 
of the rest of the posteriors. The bias on n run due to recom- 
bination uncertainties is not very significant (changes from 
n run = -0.0012 ±0.0050 to n run = -0.0034 ±0.0050, i.e. 0.4 
sigmas) but has to be taken into account. 

For completeness, we have run a 9-parameter case, in 
which we allow to vary simultaneously r, nx and n Iun in 



addition to the other six parameters. In this case, we have 
checked that the contours shown in Fig.[7]and Fig.[8]are not 
affected by the inclusion of the other parameter. 

4-1.3 Curvature 

Another possible extension of the minimal six-parameter 
model is to constrain simultaneously the spatial curvature, 
Qk- The inclusion of this additional parameter introduces 
a practical complication, since in our new parametrization 
Pncw (Eq. [2j, the Ho parameter is highly correlated with 
Q.K- Fig. |5] presents the posterior distributions for this case, 
in which the degeneration between Ho and Qk is clearly 
visible. There are two things to note. First, there is no sig- 
nificant bias to S1k due to the inclusion of recombination 
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Figure 7. Biases on the two dimensional marginalised constraints 
(68% and 95%) on inflationary parameters r-ng. Shaded contours 
represent the constraints inferred with the complete recombina- 
tion history, while the solid lines show the constraints using Rec- 
FAST vl.4.2. See text for details. 
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Figure 8. Same as Fig. [7] but in the ns-n Tun plane. 



uncertainties, as one would expect since this parameter is 
mainly constrained by information at angular scales around 
the first Doppler peak. Second, the shape of the remaining 
posteriors and the biases to the parameters are not signifi- 
cantly affected by the inclusion of this additional parameter. 



4-1-4 Residual SZ clusters/point source contributions 

One common extension of the minimal model is the inclu- 
sion of some parameters describing the residual contribution 
of Sunyaev-Zeldovich (SZ) clusters or point sources which 
are left in the maps after the component separation pro- 
cesses. For the case of Planck satellite, these are known to 
be the major con taminants at small angular scales (see e.g. 
iLeach et al]l2008h . 

For illustration, here we will consider the case of 
the residual SZ cluster contribution. One of the simplest 
parametrizations of the SZ contribution to the angular 
power spectrum is to used a fixed template from numerical 
simulations, and fit for the relative amplitude by using an 
additional parameter, -Asz- Th is approach is similar to the 
one used in lSpergel et al] (|2007l ) , where they parametr i ze the 
SZ contribution by the model of iKomatsu fc Seliakl (|2002l ) 



but allowing for a different normalization through the Asz 
parameter. 

When the parameter constraints are infered with the 
inclusion of this additional parameter, we find that (neglect- 
ing the intrinsic bias due to the degeneracy between ns and 
^4sz) the relative bias is practically not affected. In particu- 
lar, we obtain a bias of -1.42, -0.44, -2.09 and -1.03 sigmas for 
flhh 2 , Ho, ns and log(10 10 J 4s), respectively. Those numbers 
are comparable to the net biases which have been obtained 
for the minimal model in Table [2] 



5 IMPACT ON PARAMETER ESTIMATES 
USING PRESENT-DAY CMB 
EXPERIMENTS 

One would expect that the order of magnitude of the cor- 
rections to the recombination history discussed in previous 
sections (at the level of l%-2%) would have a negligible im- 
pact on the parameters constraints that we would infer from 
present-day CMB experiments. As shown in FCRW09, the 
changes introduced in the power spectra (both temperature 
and polarization) are significant at high multipoles, in the 
sense that they are larger than the benchmark level esti- 
mated as ±3,j i (see lSeliak et alj|2003h . 

To quantify this fact, we have obtained the poste- 
rior distributions for the case of the minimal model with 
six free p arameters, combining the CMB information from 
WM A.P5 dHinshaw et al1l2009h . ACBAR (iKuo et alj|2007l) 
CBI dSievers et al.ll2007f) and Boomerang ijjones et al.ll2006l : 
iMontrov et al] l200d ). together with measurements on the 
linear matter power spe ctrum based on lu minous red 
galaxies from SDSS-DR4 (jTegmark et alj|2006t ). Figure ITOl 
presents the results for the case of using the standard Rec- 
fast recombination history, together with the case of us- 
ing our most complete description of the recombination his- 
tory, as presented in the previous section. As expected, the 
modifications on the shape of the posteriors are very small 
and no biases are seen in the parameters except for ns and 
log(10 10 J 4s), which are slightly biased. Our analysis includ- 
ing the full description of the recombination history gives 
n s = 0.970 ± 0.013 and log(10 10 A s ) = 3.075 ± 0.038, while 
the result using Recfast vl.4.2 gives ns = 0.9671q 012 an d 



log(10 IU 4 s ) = 3.066 



+0.038 



In other words, this is a < 



-0.25 



and ~ —0.22 sigma bias on ns and log(10 A3), respectively. 
For completeness, we have run also the MCMC for the case 
of using WMAP5 data alone. In that case, the bias decreases 
to < —0.15 sigmas for those two parameters. 



6 DISCUSSION 

In this section, we now discuss the results presented in this 
paper focusing in three particular aspects. On one hand, we 
discuss the robustness of our results against possible mod- 
ifications of the physical description of the recombination 
process. Second, we also consider the dependence of the ob- 
tained biases if additional parameters are included in the 
MCMC analysis. Finally, we discuss the possible impact of 
recombination uncertainties on the results obtained from 
other cosmological probes different from CMB anisotropies. 
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Figure 9. Posterior distributions obtained with the remaining recombination uncertainties (see text for details) for the case of seven 
parameters, adding the curvature !2k- Solid line corresponds to the posteriors using the correct description of the recombination history, 
while the dashed lines represent the case in which the current description (Recfast vl.4.2) is used. 



6.1 Dependence of the results on the description 
of the recombination process 

As discussed above (Sect. [5]), there is a wide agreement in 
the community about the list of physical processes which 
should be included in the description of the cosmological 
recombination process. In many cases, these physical pro- 
cesses have been treated separately by at least two separate 
groups, and the agreement on the signs and amplitudes of 
the corrections is excellent in most of the cases (e.g. see the 
compilation of uncertaintie s in the physics of recombination 
in Table 2.1 of IWong||200ct ). Although an agreement at the 
level of < 0.1% is still not reached, we are almost there, as 
the remaining uncertainties seem to be at the level of 0.1%- 
0.3% between the different groups. In this sense, one would 
expect that a code which includes self-consistently all those 
processes should obtain essentially the same biases that have 
been described in Sect. [J] and have been reported in Tabled 
However, apart from the processes described in sec- 
tion 12.31 there are still some possible uncertainties which 
might lead to measurable biases on the cosmological param- 
eters. Below we now briefly address them. 



6.1.1 Hydrogen recombination 

One effect which might lead to additional biases on the cos- 
mological parameters is the inclusion of very high-n states 
in the cosmological hydrogen recombination. The computa- 



tions in this paper are based on a training set which uses 
n — 75 shells to describe the hydrogen atom. To explore the 
dependence of higher number of shells, we have repeated the 
standard six-parameter computation for the mock Planck 
data presented in Sect. 13.11 and Sect. [5] but taking as a 
reference model the one computed using n = 110 hydrogen 
shells, and trying to recover it with RlCO (which uses n — 75 
hydrogen shells). The recovered posteriors using the RlCO 
code in this case do not show any appreciable bias in any of 
six parameters of the minimal model. This is illustrated in 
Fig. 1111 where we present the posterior distribution for the 
ns parameter, which is the one having the largest bias. 

6.1.2 Helium recombination 

As mentioned in Sect. 12.31 and discussed in FCRW09, we 
do not expect major changes in our current understand- 
ing of the helium recombination, at a level which might 
be relevant for the computation of CMB anisotropics. This 
statement assumes that the uncertainties in the modelling 
of the helium atom a r e well controlled, but, as described in 
iRubino-Martm et al.l (|200et ). we still lack an accurate de- 
scription of the photoionization cross-sections, energies or 
transition rates for the Helium atom, which might lead to 
small changes in these results. One may also wonder whether 
the small differences found for the correction caused by He- 
lium feedback processes (see Sect. I2.3~2^l could matter. 

By far, the most relevant correction to the helium re- 
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Figure 10. Impact on parameter estimates for present day data (WMAP5+ACBAR+CBI+B03), together with SDSS LRGs. For this 
minimal six-parameter case, there are no obvious biases on the parameters, although there is a small difference in the shape of the 
posterior for ng. As in Fig. [6] solid lines represent the posterior distributions obtained when our most complete description of the 
recombination process is used, while dashed lines use the current description (Recfast vl.4.2). 



combination history is caused by the absorption of He i 
photons by neutral hydrogen, although the inclusion of 
the 2 3 Pi-l 1 So intercombination line also gives some con- 
tribution. These two effects are already included in the 
codes, both in R lCO (FCRW09) and in Recfast vl.4.2 
l|Wong et al.l feoolj ). For completeness, we have also quan- 
tified for this paper the impact that these two corrections 
to the Helium recombination would have on the recovered 
cosmological parameters. Our computations show that, if 
these corrections are not included (which in practise corre- 
sponds to setting RECFAST_Heswitch = in RECFAST vl.4.2, 
or equivalently, using Recfast vl.3), then the resulting bi- 
ases on the parameters are found to be -3.2, -2.0, -1.2 and 
-0.7 sigmas for ns, Qhh 2 , log(10 10 J 4s) and Ho, respectively. 
For illustration, we show in Fig. [12] a comparison of the 
posterior distribution obtained for the ns parameter when 
using the Recfast code with (RECFAST_Heswitch = 6) and 
without (RECFAST_Heswitch = 0) all the corrections to the 
helium recombination. 

This simple computation shows that additional correc- 
tions to the helium recombination history at the ~ 0.1% 
level will not matter much to the analysis of future Planck 
data. Therefore, the physics of helium recombination already 
seems to be captured at a sufficient level of precision, when 
including the acceleration caused by the hydrogen contin- 
uum opacity and the 2 3 Pi-l 1 So intercombination line, which 



together lead to a ' 
(jFendt et al]|20Qg| ). 



-3% correction to X e (z) at z ~ 1800 



6.2 Recombination uncertainties and other 
extended cosmological models 

In addition to those extended models described in Sect. 14. l"l 
there is a number of possible non-standard models for which 
the inclusion of refined recombination physics might be of 
importance. For example, when using current CMB data 
to constrain the presence of hypothetical sources of Lycv 
reso nance radiation or ionizing photons at h igh redshifts 
fe.g. lPeebles et al.ll200ol : lBean et af1l2003l . l2007h : or to probe 
dark matter models with large a nnihi l ation cross-section 
(eg. 



Finkbeinerl 120051 ; iGalli et all |2009| ; 
iHuetsi et aHl2009h ; or energy release by long-lived unstabl e 
particles (|Chen fc Kamionkowskj|20o4 IZhang et all 120071 '); 



Padmanabhan 
2009T 



or when exploring the variati o n of fundamental constants 
with time (see e.g . IGalli et alJ ll2009l) for Newton's gravita- 
tional constant, or lLandau et alJ (|200ct ) for the fine structure 
constant and the Higgs vacuum expectation value), it is ob- 
vious that neglecting physically well understood additions to 
the recombination model, as described in Sect. [5] could lead 
to spurious detections or confusion, in particular, when the 
possible effects a re already known to be rather small (e.g. see 
I Galli et alj|200g| . in the case of dark matter annihilations). 
More generally speaking, given that the largest re- 
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n s 

Figure 11. Bias on the ng parameter due to the inclusion of ad- 
ditional number of shells in the description of the hydrogen atom. 
The black solid line has been obtained from a fiducial model with 
n = 110 hydrogen shells. The dashed line corresponds to the pos- 
terior distribution recovered using the RlCO code only, with a 
training set based on n = 75 shells. The dotted line corresponds 
to the posterior distribution recovered using the current Recfast 
code. Note that we did not include for this computation the cor- 
rections due to Ly-o radiative transfer and Raman scattering. 




Figure 12. Bias on the ng parameter for different modellings 
of the Helium recombination history. The black solid line cor- 
responds to the posterior distribution recovered using the RlCO 
code and the additional corrections used in Sec. [4] The dotted line 
corresponds to the posterior distribution recovered using the cur- 
rent Recfast code (vl.4.2), while the dashed line uses the same 
version of the code but without any correction to the helium re- 
combination history. 



6.3 Recombination modelling and other 
cosmological probes 

The combination of CMB data with other datasets usu- 
ally helps to improve the parameter constraints, in some 
cases by breaking internal degeneracies which are inherent 
to the CMB data alone. One of the commonly used exter- 
nal datasets is the Bary on Acoustic Oscillations (BAOs). 
Ide Bernardis et all (|2009l ) ha ve recently shown th at a pos- 
sible delay of recombination (jPeebles et al.ll2000t ) by extra 
sources of ionizing or exciting photons leads to biases on 
the constraints from BAOs, because they largely rely on 
the determination of the size of the acoustic horizon at re- 
combination. Their conclusions can be directly translated 
here, stressing that a fully consistent combination of the 
constraints from CMB and BAOs should be done by using 
the same recombination history in both cases. 

In addition, we would like to point out that in 
principle one could reduce the uncertainty in our 
knowledge of the recombination epoch and its possible 
non-sta ndard extensions (e.g. due to an nihilating dark 
matter IPadmanabhan fc Finkbeinerl 120051 ) in two ways. 
On one hand, one could search for the imprint of the 
cosmological hydrogen recom bination lines on the CMB 
angular power spectrum ( Rubino-Martm et al.l 1200 5l : 



iHernandez-Monteaeudo et al.l 20071 ). which arises due to 
the resonant scatter ing of CMB p hotons by hydrogen 
atoms at each epoch (jBasu et al.ll2004l ). On the other hand, 
one could also try to directly observe the photons that 
are emitted during the recombination epoch. Today these 
photons should still be visible as small distortion of the 



bands ( 


e.g. see Dubrovichl 


19751; Dubrovich & Stolvarov 


19971: 


Rubino-Martm et al. 


20061: Chluba & Su 


lyacv 


2006al; 


Chluba et al. 


20071; iRubino-Martm et all 


20081; 


Chluba & Sunvaev 2009bl). 


These observations not only 



would open another possibility to determine some of 
the key cosmological parameters, such as the primordial 
helium abundance, the number density of baryons and 
the CMB monopo l e temp erature at recombination (e.g. 
I Chluba fc Sunvaevl l2008ai ). but they would also allow 
us to directly check our understand of the recombina- 
tion process and possib le non-standard aspects (e.g. see 
ISunvaev fc Chlubal [20091 . for an overvie w), for example, in 
connec tion with early energy releas e (IChluba fc Sunvaevl 
l2009dl ). or dark matter annihilations (lChluball2009f T 



combination uncertainties are obtained for ns, Qhh 2 and 
log(10 10 ylg), one can say that any additional parameter 
showing a strong correlation with those three might be 
biased if an incomplete description of the recombination 
physics is used. In this sense, neglecting the refinements 
to the recombination model could be as important as not 
taking into account, for instance, uncertainties in the beam 
shapes, which also have been shown to co mpromise our abil- 
ity to measure ns (|Colombo et al.ll2009T l; or the combined 
effect of beam and calibration uncerta inties, which intro - 
duce significant biases to ns and Q,\,h 2 l|Bridle et al.ll2002T ). 
although other parameters (like Qk) are essentially not af- 
fected because these are basically constrained by the posi- 
tion of the peaks, and not by their amplitudes. 



7 CONCLUSIONS 

In this paper, we have performed a MCMC analysis of the 
expected biases on the cosmological constraints to be derived 
from the upcoming Planck data, in the light of recent de- 
velopments in the description of the standard cosmological 
recombination process. Our main conclusions are: 

• An incomplete description of the cosmological recom- 
bination process leads to significant biases (of several sig- 
mas) in some of the basic parameters to be constrained by 
Planck satellite (see Table[5}, and in general, by any future 
CMB experiment. However, these corrections have a minor 
impact for present-day CMB experiments; for instance, us- 
ing WMAP5 data plus other cosmological datasets, we find 
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a ~ —0.25 and ~ —0.22 sigma bias on ns and log(10 10 As), 
respectively, while the rest of the parameters remain un- 
changed. 

• Today, it seems that our understanding of cosmological 
recombination has reached the sub-percent level in X e at 
redshifts 500 < z < 1600. However, it will be important to 
cross- validate all of the considered corrections in a detailed 
code comparison, which currently is under discussion among 
the different groups. 

• Given the range of variation of the relevant cosmological 
parameters, it is possible to incorporate all the new recombi- 
nation corrections by using (cosmology independent) fudge 
functions. Here we described one possibility which uses a 
simple correction factor to the results obtained with Rec- 
fast (see Sect .[3~2")l . We provide the function f(z) on the 
RlCO-webpage Q 

• The physics of helium recombination already seems to 
be captured at a sufficient level of precision, when including 
the acceleration caused by the hydrogen continuum opacity 
and the 2 Pi-1 So intercombination line. The biases caused 
by neglecting only these corrections are -0.8 and -0.4 sigmas, 
for ns and Qbh 2 , respectively. 

• When allowing for more non-standard additions to the 
recombination model (e.g. related to annihilating dark mat- 
ter), the biases introduced by an inaccurate recombination 
model could lead to spurious detections or additional con- 
fusion (see Sect. 16. '2\ . 
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